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DETERMINATION OF GENERAL RELATIONS FOR THE BEHAVIOR OF TURBULENT 

BOUNDARY LAYERS 

By Albbbt E.-von Dobnhoff and Nbai> Tbteevin 



SDMMAEY 

An analysis has been made of a considetable amount oj data 
for iurhvleni bowndary layers along wings and bodies of various 
shapes in order to determine the fwndamenial variables that con- 
trol the development of turbulent boundary layers. It vxls found 
that the type of velocity distribution in the boumdary layer could 
he expressed in terms of a single parameter. This parameter 
was chosen as the ratio of the displacement thickness to the 
momentum thickness of the boundary layer. The variables thai 
control the development, of the turbulent boundary layer appar- 
ently are (1) the ratio of the nondimen-sional pressure gradient, 
expressed in terms of the local dynamic pressure outside the 
boundary layer ajid boundary-layer thickness, to the local skin- 
friction coefficient and (5) the shape of the boundary layer. An 
empirical equation has been developed in terms of these variables 
that, when used with the momentum equation and the skin- 
friction relation, makes it possible to trace the development of 
the turbulent boundary layer to the separaHon point. 

INTRODUCTION 

A good measure of the understanding of the general prob- 
lem of the flow of a real fluid about a body of arbitrary shape 
can be taken as the degree of approximation with which the 
aerodynamic characteristics of tiie body can be calculated 
on the basis of existing knowledge. The flow in regions re- 
moved from the wake and from the surface of a body obeys 
very closely the laws for perfect frictionless fluids at all rea- 
sonably large Reynolds numbers. This type of flow is well 
understood, although the detailed computations may be 
difficult in some cases. Departure of the flow of real fluids 
from that of the frictionless fluid is caused almost entirely 
by the failure of the idealized fluid to reproduce the actual 
flow conditions at the surface of the body. 

Because the theory of perfect fluids gives zero drag for all 
bodies and gives no information concerning conditions that 
lead to separation of the flow from the surface, these pheno- 
mena must be almost entirely associated with the behavior of 
the flow at the surface or, in other words, with the boundary 
layer; that is, all drag, with the exception of induced drag 
and drag due to shock waves, all cases of flow separation such 
as occur at maximum lift or at high aileron deflections, and 
all Reynolds number effects are entirely dependent on the 
behavior of the boimdary layer. 

Because the flow-of ideal fluids is well tmderstood, the 
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problem of calculating the actual flow resolves itself into a 
study of boimdary layers. There are three general t3rpes of 
boimdary layer: laminar, transitional, and turbulent. In 
most cases, the important laminar-boundary-layer character- 
istics — thidkness, skin friction, and point of separation — can 
be estimated with sufficient accuracy from methods described 
in references 1 and 2. If , in unusual cases, more detailed in- 
formation is required, one of the methods of calculation given 
in reference 3 may be used. 

Experiments in low-turbulence wind tunnels and in flight 
have shown that extensive laminar layers may be maintained 
in the presence of a favorable pressm-e gradient. Although 
little is known about methods of determining the position of 
the transition point, transition must occur either associated 
with laminar separation or at some position upstream of the 
separation point. The position of the transition point varies 
widely within this region because of changes in stream tur- 
bulence, surface condition, and Reynolds number. 

A considerable amoimt of progress has been made in de- 
termining turbulent sMn friction in pipes and along flat 
plates with zero pressure gradient. As a result of the de- 
velopment of the "mixing-length" theory, notably by 
Prandtl and von Efiimdn, reHable rules that should be ap- 
pHcable throughout an extremely wide range of Reynolds 
numbers have been found for calculating the skin friction 
and velocity distribution in the boundary layer along flat 
plates and in pipes. 

In order to flnd the effect of varying pressiu:es on the 
boundary-layer characteristics, von E!£rm&n applied the 
momentmn theorem to the boundary layer and derived the 
so-caUed "momentum equation" that gives the rate of 
thickening of the boundary layer if the type of velocity dis- 
tribution within the boundary layer, the external pressure 
distribution, and the skin friction are known. It was found 
in a number of cases that, if the shape of the boundary layer 
and the sMn friction were assumed to be the same as in 
pipes, good agreement was obtained between the calculated 
and experimental boundary-layer thickness and skin friction. 
This procedure, however, faUs to give any information con- 
cerning the changes in boundary-layer characteristics that 
lead to separation of the flow from the surface. In addition 
to the momentum equation, a relation is needed between the 
shape of the velocity distribution in the boundary layer, the 
skin friction, and the pressure distribution. 
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Several attempts have been made to find such, a relation. 
From experiments in converging and diverging channels, 
Nikuradse, using water (reference 4), and Donch, using air 
(reference 6), fotmd that, in cases in which the boundary 
layer along the walls met in the center of the channel, the 
velocity distribution across the channel was a fimction of 
ari/R, where a is the angle of divergence of the channel and 
jB is the Reynolds niunber based on the channel width and 
average velocity. Buri, whose work was discussed by 
Prandtl in reference 6, used the results of Nikuradse and 
DSnch in an attempt to calculate the general behavior of 
turbulent boundary layers.- Buri assumed that the shape of 
the boundary layer was always given by the parameter of 
Nikuradse and Donch. One of the weaknesses in Buri's 
calculations was the assumption that the shape of the 
velocity distribution depended only on , the value of this 
parameter and was independent of the previous history of 
the boimdaxy layer. It may be pointed out that the para- 
meter of Nikvuradse and Donch is not essentially a shape 
parameter but represents a function of the pressiure gradient 
and skin friction which is assumed to determine the shape of 
the velocity distribution. 

Gruschwitz (reference 7), using a parameter dependent 
upon the type of velocity distribution in the boundary 
layer, found a relation between the pressure gradient and 
this parameter. This relation, with the momentum equa- 
tion, was sufficient to determine the development of the 
turbulent boundary layer along a surface. Although 
Gruschwitz obtained good agreement with experiment for 
the data presented in his paper, other investigators whq have 
tried to use the method reported poor agreement forcasesin 
which the turbxilent boundary layer separated from the sur- 
face. Peters (reference 8) conducted an investigation for 
the specific purpose of testing the Gruschwitz method of 
calculation and concluded that the Gruschwitz method 
cannot be used to determine the location of the separation 
point nor even, in many cases, to predict whether separation 
will occur at aU. 

The purpose of the present investigation is to determine 
the important variables that control the behavior of turbu- 
lent boimdaiy layers and to develop general relations in 
terms of these variables that describe the boundary-layer 
motion. The fundamental variables must, of course, be 
expressed nondimensionaUy in terms of local boundary-layer 
quantities. 

The quantities at a given station along a surface that were 
felt to have the most important effect on the further develop- 
ment of the boundaiy layer are the following: 

(1) Shape of the boundary-layer profile 

(2) Rate of change along the smJace of the dynamic 

pressiu:e outside the boundary layer 

(3) Skin friction 

Data from various published somces and from tests in the 
NAOA two-dimensional low-turbulence timnel made specifi- 
cally for the present investigation were analyzed in terms 
of the foregoing variables in order to find the needed general 
relation for the rate of change of shape of the boundary- 
layer profiles. 
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SYMBOLS 

angle of divergence of channel; also, angle of attack 

Reynolds number 

minor axis of ellipse 

velocity within boundary layer 

free-stream velocity 

distance along surface 

chord 

distance perpendicular to surface 

velocity outside boundary layer 

skin friction per unit length 

dynamic pressure outside boundary layer 

shape- parameter 

displacement thickness |^J^ ^1 <^2/J 
density 



momentum thickness 
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"Viscosity 

boundary-layer thickness 
kinematdc viscosity- 
constant 

equivalent length of flat plate before pressure recovery 

begins 
width of channel 

total pressure at distance 6 from surface 

free-stream total pressure 

shape parameter j^l ~(^^^ J 

-velocity at distance 6 from surface 

initial value of g 

constant 

initial value of 0 

initial value of Be 



5^=2.557 log. 4.075 Be^ 
^■=2.557(^-^1) 



EXPERIMENTAL DATA 

The data used in the analysis were collected from the 
available literature and from tests performed in the NAOA 
two-dimensional low-turbulence tunnel. The following table 
sho-wB the data used in the present analysis: 
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10.1 


Present report 
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Present report 
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The paper by Gruschwitz (reference 7) contained data on a 
Vidng and on channel walls obtained from tests at Gottingen. 
The data from run 2 on the channel wall in the form of 
boundary-layer velocity profiles, pressiu-e distribution, and 
curves of the momentum thickness and the shape parameter 
plotted against distance along the plate were used. From 
Peters' paper (reference 8), data were taken in the form of 
curves of the momentum thickness, shape parameter, and 
pressure coefficient against the position along the airfoil 
chord. These tests were made in order to check the method 
of calculation proposed by Gruschwitz. The [Reynolds 
number of the test was not given explicitly and was judged 
to be slightly less than 4,000,000. The results from the tests. 
at an angle of attack of 9° were used because the data for 
this angle of attack were presented in a convenient form and 
because separation of the flow had taken place at the rear of 
the wing. 

A few points were obtained from the boundary-layer veloc- 
ity profiles and pressure distribution contained in reference 11. 
The data were obtained from a test of an elliptic cylinder 
afc an angle of attack of 0° and a Reynolds mmiber of 
118,000 based on the minor axis D of the ellipse. The ratio 
of major to minor axis was 2.96, making the Reynolds num- 
ber based on the major axis equal to 349,000. 

The data on turbulent -boundary layers involving separa- 
tion obtained from the NAOA two-dimensional low-tiurbu- 
lence tunnel were collected from previous tests and from 
tests performed specifically for the present investigation. 

A few turbulent-boiindary-layer profiles that were. not 
close to separation were obtained from previously published 
data on the NACA 0012 airfoil (reference 10). These points 
were used in the analysis mainly beC|ause the boundary layer 
was far from separation and the use of these points helped 
to give a better distribution of data. The data obtained 
on the NACA 66,2-216 airfoil are given in reference 9. 
Ordinates for this section can be foimd by methods described 
in reference 12. 

Tests of a thick airfoil, the NACA 65(216)-222 (approx.) 
were made at three Reynolds numbers at two angles of 
attack. The methods of obtaining the data were the same as 
those described in reference 9. All the tests of this airfoil 
involved tiu:bulent separation. The data from these tests 
are presented in the form of boimdary-layer velocity profiles 
for a number of stations along the chord (figs. 1 to 4). The 
pressure distributions are given in figiu-e 5. The region of 
turbulent separation is indicated in the pressure distribution 
as the flat region at the rear of the airfoil. The chord of the 
airfoil was 24 inches and the airfoil was finished as described 
in reference 9. The ordinates can be derived by methods 
explained in reference 12. The finiah was free from all 
surface imperfections that could be felt by hand but had a 
strip of carbonmdiun-covered cellulose "Scotch" tape 1 inch 
wide on the upper surface near the leading edge. 




(8) B-.O.MX10». 
JJ-1.61X10«. 
I JJ-2.67X10». 

FlQtlKKl.— Bonndary-layer velodty profiles In the resion forward of the 36-peiccntHiord 
station. Airfonseoilon, NAOA 65(216^^222 (approx.); acg.!". 
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(a) £>=>0^X10<. 

(b) iJ-lJlX10«. 
(o) jB-2.67X10«. 

FiQTTBE 2.— Boimdaiy-IaTei Telocity profiles in the region Irom the 35-peioent-cbord to the 65-percentohord stations. Airfoil section, NAOA e6(2ie)-223 (approx,); a-8.1°. 
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The other model tested in the NAGA two-dimensional 
low-turbulence tunnel for data to be used in the present 
analysis was the NACA nose-opening airfoil shape 13. The 
model had a chord of 36 inches and was finished iu the same 
manner as the NACA 66,2-216 (reference 9). The section 
ordinates for the NACA nose-opening airfoil shape 13 are 
given in reference 13. The wing was tested at an angle of 
attack of 9.1° at three Ke3Tiolds numbers. The turbulent 
separation obtained in this test was not so marked as that 
obtained in the tests of the NACA 66,2-216 and NACA 
65(216)-222 (approx.) airfoils, although tufts placed at the 
rear of the wing on the upper surface indicated separation. 
The data from these tests are presented in figures 6 and 7 in 
the form of boundary-layer velocity profiles for a mmiber of 
stations along the chord. The pressure distributions are 
given in figure 8. The beginning of separation is indicated 
in the pressure distributions by the flattening of the curves 
at the rear of the airfoil. The small flat region in the pressure 
distribution at the nose of the airfoil is an indication of lami- 
nar separation. The boundary-layer velocity profiles for 
the region at the nose of the airfoil are^own-in'-figiure 6. 
The peculiar shape of the velocity distribution for some of 
the stations, particularly in the curves that shew^increasing 
velocity ^vith approach toward the wall, is probably caused 
by spanwise flows over the airfoil. The boundary-layer 
thiclmesses obtained in these tests were much larger than 
those usually obtained for airfoils of 36-inch chord. 

ANALYSIS 

The equation that gives the rate of change of the momen- 
tum defect in a boimdary layer, originally derived by von 
KArmdn, may be written in the following form for two- 
dimensional flow: 



dJd ■ / g+2 \ edq To 



whore 



d momentum thickness j^J" -^^i— dyj 

u velocity within boundary layer 

U velocity outside boundaiy layer 

y distance perpendicular to surface 

To sMn friction per luiit length 

g dynamic pressure outside boundary layer 

X distance along surface 

II shape parameter 

The shape parameter H is defined as the ratio 6*/tf where 

5*= (^—^ dy. The difference between the actual flow 

of momentum in the boundary layer and that of the same 
quantity of fluid flowing with velocity U is plPd. From this 
relation, the length d is given the name momentum thickness. 




0. 53X10<. 

1. EIXIO*. 
3.MX10*. 

FiointE 3.— Boundary-layer Tdadty profiles in the region forward of the SE-i>erceiit-cli<iTd 
utatlon. AJrfoD section, NACA 66(21(0-222 (approx.); a-I0.1°. 
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(a) ii-0.93X10i<. 

(b) JJ=1.HX10». 
(o) ie-2.MX10». 

FiQVBB 4.— Boondaiy-Iayer velocity proffles fa the region from the 35-pcrcent-chord to the dS-percent-chord itatlons. Alifoll section, NACA 66(216)-222 (appiox.); aolO.l 
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FioCB^ S^XTppor-sortace pressure distrtbutions for several BeynoMs nombeia. Note- 
opeoini; airfoil sbape 13; a-9.1<>. 

The length 5*, called the displacement thickness, is the 
amount by which the streamlines just outside the boimdary 
layer are displaced because of the reduction of velocity 
within the boundary layer. Because 6 depends on the second 
power of the velocity distribution, whereas 5* depends on 
only the first power, the ratio h*lO depends on the maimer in 
which ujU varies with y — that is, upon the shape of the 
boundary-layer profile. 

The momentxim equation in the form just given contains 
only local boundary-layer quantities. The local skin-friction 
coefficient is To/g. The nondimensional pressiire gradient is 

^ ^> where 6 is the unit of length and g is the unit of dynamic 

pressure. 

Although it has been shown that the. shape of the bound- 
ary-layer profile determines M, the converse cannot be 
proved from mathematical considerations alone. If H does 
actually determine the shape of the boimdary-layer profile, 
then all points of ujU plotted against H at a constant value 
of yjd should fall on a single curve. A collection of such 
curves for various values of ylQ is shown in figure 9. The 
data presented in figure 9 represent the collection of all the 
boundary-layer profiles that enter into the analysis. Figure 9 
shows that ujU is a function of M alone for a given value of 
yjd. This conclusion is important because it means that 
turbulent-boundary-layer profiles form a single-parameter 
family of curves. The complete velocity distribution in the 
boundary layer is known when 6 and Hhave been determined. 
Figure 10, which is a crpss plot of figure 9, gives turbulent- 
boundary-layer velocity profiles corresponding to Various 



values of ff. As the separation point is approached, the 
value of M increases. Because the tmrbulent separation point 
usually is not very weU defined, it is not possible to give an 
exact value of H corresponding to separation. The value of 
H usually varies so rapidly near the separation point, how- 
ever, that it is not necessary to fix accuratdy the value of H 
corresponding to separation. Separation has not been 
observed for a value of if less than 1.8 and appears definitely 
'to have occurred for a value of M of 2.6. Gruschwitz's 
criterion for imminent separation is equivalent to a value of 
M of 1.85. 

The fact that the type of velocity distribution in the 
boundary. layer can be given in terms of a single parameter 
greatly simplifies the study of turbulent boundary layers. 
It is now necessary to determine only the manner in which 
this parameter varies along the surface as a function of the 
external forces acting on the boundary layer. 

The external forces acting on the boundary layer at any 
point are the pressure gradient, expressed nondimensionally 

as and the skin friction, expressed nondimensionally 

as To/ff. The assumption is made that the rate of change of 
H rather than H itself is related to the local forces. This 
assumption is dearable in order that the boundary-layer 
conditions downstream from a point shall be definitely con- 
nected with the conditions upstream of the point; that is, a 
sudden change in the pressure gradient should not produce 
a discontmuity in the type of velocity distribution in the 
boundary layer. Prandtl (reference 6) has pointed out this 
difficulty both in the Fohlhausen theory for laminar bound- 
ary layers and in Buri's method for calculating turbulent 
boimdary layers. Expressed in nondimensional form, the 

rate of change of H is given as B 

In tiie early stage of the analysis, the experimental data 

c[dx 

tion was obtained for a limited amoimt of data. As the 
analysis was extended to include more data, systematic 

variations with Reynolds number were noticed. When 0 

gdXTo 

Reynolds number was eliminated. The skin friction was 
tentatively assumed to be given by the Squire and Toung 
formula (reference 14) 



were plotted in the form of 6 against - • 



was plotted against -r ^ the consistent variation with 



where 



^=[5.890 logio (4.075S«)]* 
To 



This formula was chosen because of the good agreement 
obtained between the experimental drag coefficients for air- 
foils and those calculated by the Squire and Young method. 
It was felt that the local skin-friction coefficients thus 
determined were probably more accurate in most cases than 
those determined directly from the l)oundary-layer surveys. 
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The quantity ^ ^ qj. a quantity differing only by a 

constant factor, has frequently been used as a parameter 
for boundary-layer phenomena. • For example, the Pohl- 

hauson parameter for laminar boundary layers ^ where 
5 is the boundary-layer thickness SLadj) the kinematic viscos- 
ity, can be shown to be equivalent to ^ For a lami- 
nar boundary layer, to=^ (^) (^) Proportional to 
U/5. It thus follows that /i/5 is proportional to to/Z7. Sub- 
stituting for in the Pohlhausen parameter and replao- 
pU po 

ing ^ by its equivalent ^ ^ gives a quantity proportional 
to the Pohlhausen parameter O -1^^, which ia turn is 



By a similar process 



equal to some constant times ^ 

of reasoning, the parameter where « is the equivalent 

length of the flat plate before pressure recovery is begun, 
which determines the amount of pressure that can be re- 
covered in a laminar layer with a straight-line velocity 
gradient as given in reference 1, can be shown to be propor- 



tional to — ^ 
2 ox To 



Nikuradse (reference 4) found that his re- 



sults for turbulent flow, which ^ve the velocity distribution 
across divei^ing or converging channels, agreed with similar 
measurements made by Donch (reference 5) when a had 
the same value. For a given type of velodty distribution 

across the channel, ^ ^ *t the center is proportional to the 

angle of divergence, where b is the width of the channel. 
"Within the range of Reynolds numbers covered by the 
investigations of Nikuradse and D5nch, the skin-friction 
coefficient ro/g at the wall was inversely proportional to -v^- 
For a given type of velocity distribution, a is therefore 

proportional to ^ ^ It may be pointed out that, 

although in Buri's theory of tvurbulent boundary layers and 



Pohlhausen's theory of lanunar boundary layers 



was 



g_dx Tfi 

assumed to determine the type of velocity distribution in the 
boundary layer, in the present analysis it is assumed that 

- ^ ^ affects only the rate of change of the type of velocity 
2 ox To 

distribution. 

It seemed highly probable that the rate of change of H 



should depend not only on the ratio of the pressure gradient 
to the skin friction but also on the value of H itself. Plots 

were therefore made of against -4^^ at 0.1 intervals 

ax 2 ax To 

of H for aU the data entering into the present analysis. 

These plots are given in figure 11. Although the points 

show considerable scatter, definite trends for the variation 

of with both - ^ ^ and H are observable. It may bo 
ax 2 cte To 

pointed out, however, that both ^ and ^ were the slopes 

of experimentally determined curves. Large scatter of the 
data therefore is to be expected. The laige scatter of the 
points in figure 11, consequently, does not necessarily in- 
dicate any serious inadequacies in the present analysis. 
From a study of the available data, it was found that the 



. ^ , JH ... edq2q 
vanation of O-j- with - ^ 
ax 2 oa; To 

represented by the equation 



and H cotild bo fairly well 



eg=g4.680(H-»..7« [-^ ^ ^-2.036(£r-1.286)] 

The exponential form of the factor multiplying the second 
member of this equation was chosen because the data for 
high values of H, although not very complete, nevertheless 

indicated that ^ was large. It may be noted that 1.280 



dx 

is the value of H for the =-power distribution of velocity in 



the boundary layer. It is seen from the foregoing equation that 
no change in Hia indicated for the case in which ^=0 and H 
has the value 1.286. The degree to which the equation for 
approximates the experimental data from which it was 
derived may be seen in figmre 11. Each of the straight linos 
in figure 11 was obtained from the equation for by giving 

H successive values of 1.35, 1.45, 1.55, and 1.65. The slope 
of these lines is given by the factor e*-"^*^"*-*^", and the inter- 
cept for 0^=0 is given by the term —2.035 (S— 1.286). 



COMPARISON OP PRESENT ANALYSIS WITH GRUSCHWTZ 

ANALYSIS 

In reference 7, Gruschwitz analyzes the behavior of turbu- 
lent boundary layers in the following manner: 



2 



'Fiv, Be) 



where 

1} shape parameter j^l i"^ the same sense as H 

ue velocity at y=0 

and 

hi total pressiure at distance 6 from surface (gi in reference 7) 
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(c) Valnesof i7iTom Ui to 1.6. 
(Si Vataes otHtatm 1.6 to 1.7. 

FiQDBE 11. — Oonclnded. 
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In writing this relation, it is implicitly assumed that 6 ^ is 

independent of Because of this assumption, the Gru- 

schwitz analysis is subject to the same critidsm as are the 
Buri and Pohlhausen methods for making boundaiy-layer 
calculations. Now 

r]g=ho — hi 

•where Ao ((/o in reference 7) is the free-stream total pressure; 
hence, 

dx " ''"^ 



2 dx 



2 dx 



and therefore 



Except for the factor 2g/To, this relation is similar in form 
to the relation found from the present investigataon; namely, 



The Gruschwitz relation, however, is very restricted in the 
form of the dependence ^ ^ ^ comparison with 
the type of relation used in the present investigation. Fur- 
thermore, no variation of fl^^ indicated in the 

final equations given by Gruschwitz. 

In the Gruscmvitz analysis, the arbitrary function wliicli 
was to be determined experimentally contained only one 
variable ^; whereas, ia the present analysis, the arbitrary 

function contains two variables, - ^ ?2 Apart from 

' 2 (KT To 

the neglect of Be, one reason for the failure of the Gruschwitz 
analysis is that a correlation of all tuxbulent-boundary-layer 
data in terms of a function containing only one independent 
variable was not possible. 

METHOD OP CALCULATION 

For calculating the characteristics of turbulent boundary 
layers, the following information is required: The initial 
values of 6 and H, the pressure distribution oyer the body, 
and the Reynolds number. The equations that are used in 
making a computation are 

dB H+2 6 dg^To^ - 
da;"*" 2 ? die 22 

eg=c*-'»«'^-"'« [—^^^-2-035 (ff-1.286)] 

^=[5.890 logio (4.075 5,)]* 

In order to reduce the work of computation, to/22 was plotted 
against Be in figure 12 and the factor e*-9«K»-s-wo ia the 

dS 

equation for was plotted against H in figure 13. 

The momentum equation and the equation for ^ are 
simultaneous fitrst-order differential equations that can be 

740023 — 48 2(i 



solved by a step-by-step calciilation. It is usually necessary 
to use sucb a method although, for some particular cases, 
the equations may be integrated directiy. The method of 
calculation is as follows: The values of the variables entering 
into the computation at the initial station are substituted 

dS 

in the momentum equation and the equation for Values 

for ^ and are thus obtained at the initial station. An 

increment of the length along the surface of the body x is 

then chosen and multiplied by ^ and to give M and AH, 

respectively. These increments of B and H are added to the 
initial values and result in values of 0 and M for the new 
value of X. The process is repeated until the desired result 
has been attained. Separation may be considered to have 
occurred when H rises to about 2i6. 

The choice of the increment of x is a matter for the judg- 
ment of the individual investigator. As a general rule, the 

increments of x should be made small when ^ or changes 

rapidly from one value of a; to the next. In order to decrease 
the length of the calculation, the increments of x must be 
chosen as large as is compatible witb the accuracy desired. 
For the computations that were made in order to check the 
method of calculation, the increment of x for one step was 

so chosen that Aa;'<0.0025, where A ^^is the change 

in between two successive values of x and Aa; is the incre- 
ment of X. This criterion furnishes a measiuro of the maxi- 
mum error that can be expected in AH" for one step of the 
computation. When the flow approaches separation, H 
usually increases very rapidly and, in such cases, the fore- 
going criterion may be disregarded without appreciable error 
in the position of the separation point. By disregarding the 
criterion when the flow is close to separation, the length of 
the computation may be reduced. 

A sample calculation for the NACA 65(216)-222 (appros.) 
airfoil section at a=10.1° and at iJ=2.64X10° is given in 
table I. 

If the question of separation is not involved and if the 
variation of H along the siuiace is not of interest, reasonably 
accurate values of d may be obtained by assuming a constant 
value of H and merely using the momentum equation to- 
gether with the skin-friction relation to determine 6. This 
procedm-e is substantially the same as that of reference 14 
where a constant value of J? of 1.4 was chosen for calculating 
the profile drag of airfoil sections. 

TESTS OF METHOD OP CALCULATION 

In order to obtain a general check of the method of calcu- 
lation and to determine whether the scatter of the points in 
figure 11 was primarily due to the difiBculty in obtaining the 
slopes of experimental curves or to serious inadequacies in 
the analysis, computations were carried through for eight 
cases. For aU these computations, the initial values of H 
and 0 were obtained from experimental data. 
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(a) Sliapo iiammoter 

(b) MoTDOTitam thlcJcncsa fl/e. 

FlQUBE 14.— Variation of experimental and ca1mlRt4><i valaes of tbape iwranifltfr H and 
momentom tiilckness «/e iritli x/e. NAOA 66^216, iv=^.% airfoil section; Ji-2.6X10S; 

A comparison between the calculated and ezperimeatal 
variations of H and Q along the surface is shoTSTi in figures 14 
to 24. A computation made by the Gruschwitz method 
(from reference 8) is included in figure 17. For some of the 
cases, comparisons of the calculated and experimental 
boundary-layer profiles at one or two positions are also pre- 
sented. In general, the calculations are in good agreement 
with the experimental curves. No systematic differences 
were foimd between the calculated and experimental curves 
of B.. Although the agreement between the calculated and 
experimental curves of Q (figs. 14 to 24) is good in most cases, 
some consistent diOterences are apparent as the separation 
point is approached. In this r^on, many of the calculated 
values of Q are less than the experimental values. One ex- 
planation of the discrepancy, of course, is that the Squire 
and Young sMn-friction relation is in error in not indicating 
an increase in skin friction as the boundary-layer velocity 
profile approaches the shape for separation. This tendency, 
however, is contrary to the general impression that the skin 
friction should decrease as the separation point is approached. 
As the flow approaches separation, the fluctuations in the 
direction of flow increase. Such fluctuations make a pitot 
tube read velocities higher than the actual velocities. These 
fluctuations are a large proportion of the mean flow dose to 
the surface where reversed flow first begins. This behavior 
of a pitot tube may explain why the turbulent velocity pro- 



files, which are close to the separated state, all have the 
characteristic hump at small values of y/c. The velocity 
profiles for large values of H are therefore in error in the re- 
gion dose to the surface. The error in profile shape affects & 
. as wdl as ff. Reading vdocities too high for all points where 

■^^0.5 makes the integral for d, which is | 

larger than it actually is. It is to be expected, there- 
fore, that values of d determined by pitot-tube readings 
should be higher than the true values under conditions of 
unsteady flow. On the other hand, as separation is ap- 
proached, the relativdy greater velocity fluctuations near 
the surface may cause the skin friction to be higher than 
when conditions are far from separation ; and the effect of the 
lower average vdocities near the surface, such as occm* for 
higher values of H, may thus be compensated and possibly 
overbalanced. An attempt was made to correlate the local 
skin-friction coeffidentwithff, but no consistent results were 
obtained. Although there still is considerable doubt concern- 
ing the true value of the skin-friction coefficient for con- 
ditions approaching separation, it is interesting to note that 
the Squire and Yoimg skin-friction relation was used through 
a range of Rb from 600 to 48,000 and apparently gave good 
results for most of the region covered by the turbulent 
boundary layer. 




(a) Sliapo parameter JET. 

(b) Momentnm thickness 6t centimeters. 

FiaimE IS.— Variation of experimental and calcnlated valaes of sliapo parameter H and 
momentmn tliiclme^a B with x tac Sim in a nlmnnel. (Data lirom reference 7, ran 2.) 
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KlQUBE IS,— Comparlsoa of experlmontal and calcalatcd volodtr piofUes. (Data lirom 

reference 7, ran 2.) 
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(a) Shape parameter H. 

(b) Momentom thickness 

FiauRE 17.— Variation of experimental and calcniated Talnes of shape parameter H and 
momentmn th'i'ln'PCT «/e vrlth x/c for symmetrical alrfofl. R^ZXKSIIfi; a— 0°. (Data 
from roferenoe 8.) 
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FiooBE I8j— Variation of experimental and calcolated -valoes of shape parameter H and 
momentnm tlilckness with xtD tac elUptio cylinder. Jiaii6,000, based on minor axis. 
(Data from reference 11.) 



1.0 



.8 



.4 



.2 



0 
■1.0 

.8 



.6 













































































































































































































inn 
Jot 


Ten 

°d 






























— L 

— C 


al- 










i 


































































































































































































































































































































— t 


'xpt 


irin 
j/at 


Ten, 
3d 


ol- 




























-C 












/ 


t 


































t 


































.)- 





































.0/ .02 .03 .04 .05 .06 .07 X)a J09 

(a) 2-2-"- 

(b) f-X88. 

FiOTJEE 19. — Oomparlson Of experimental and calcolated velocity profiles fbr the efliptio 
cylinder. R^lViffA, based on minor axis. (Data from reference II.) 
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FlQDBE 20.— Variation of experimental and calculated valaes of shape paiametcr H and 
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FIOVBE 21. — Comparison of experimental and calculated veladty profiles for iioseK>peiiing 



In making a calculation, the initial condition of tiie bound- 
aiy layer must be known. Computations whieh have been 
made do not indicate that the calculation for the H-curve is 
espedaUy sensitive to the initial value of B. When the calcu- 
lation is to be made for a case in which the boimdary layer 
is in a strong adverse pressure gradient — that is, when 

- ^ ^ is of the same order or greater tiian2.035(fl'— 1.286) — 

2 (to To 

the initial value of H must be accurately determined. This 
fact is obvious because the value of H determines how dose 
the flow "is to separation. If the calculation is begun in a 

region where ^ ^ very small or positive, the boundary 
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(a) Shape parameter H. 

(b) Momentum thlofcoess #/c. 

TiooKE 21— Variation of experimental and calculated values of sliapo parameter // and 
momentum thickness eje wIS x/e. Airfoil section, NACA e5(3ie)-232 (opprox.)i i2-0.92X 
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FiauBE 23— Variation of experimental and calculated values of shape parameter // and 
momentum thi<dcness0/e'«rItha/& Airfoil section, NAOA 6S(316}-222 (approx.): iJ-2.67X 
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layer is not very sensitive to the initial value^of H. For 

example, if ^=0, H" will eventually have the value of 1.286 

regardless of its initial value. 

The boundary-layer thickness is not particularly sensitive 
to the initial value of H. This effect is easily vraified from 
the form of the momentum equation. The effect of changes 
in the initial value of 6 on the boundary-layer thickness 
farther dovmstream depends on the relative magnitude of 

^ and If the pressure gradient is large in comparison 

with the skin friction, a change in the initial value of 6 will 
produce a proportional change in the subsequent values of 0; 
whereas, if the pressure gradient is small in comparison with 
the skin friction, a change in the initial value of d will produce 
a constant increment equal to the initial change. 

The initial value of 6 ordinarily may be taken the same as 
the value for the laminar boundary layer at the transition 
point. Not enough is known about the mechanism of tran- 
sition to be able to state in general what the initial values 
of H should be. If transition occurs in a zero or "favorable" 
pressure gradient or if the boundary layer is sufficientiy 
thin at the transition point, in accordance with the foregoing 



discussion, the position of the turbulent separation point 
will not be greatiy affected by the choice of the initial value 
of B". 

PRESSnSE RECOVERY WITHOUT CHANGE IN BOUNDARY- 
LAYER SHAPE 

The equation for indicates that, for each value of H, 

■ 6 do dJS. 

such a value may be assigned to — ^ that d ^= 0. Pressure 

may be recovered, therefore, without a change in boundary- 
layer shape if the correct pressure distribution is used. The 
necessary pressure distribution can be obtained by using 

the equations for and ^* These equations can be 

int^ated directiy if jff is assumed to have a constant value. 

Wheng=0, 



■il^=-2-035(fi^-l-286)=C7, 



(1) 



where (7i is a constant. The momentum equation then 
reduces to the form 

where 

^7=1-1-2.035 (^^) (H-1.286) 

Elimination of T<i/2g between equations (1) and (2) gives~the 
relation 

Upon integration, 

where 6o is the initial value of d and go is the initial value of 
2. The relations between 0 and x and, consequently, be- 
tween 2 and x are obtained from integration of equation (2). 

For To/22, the Squire and Young skin-friction relation will 
be used 

1 

22 [2.657 log. (4.075 R,)]' 



where Bsf, is the initial value of B$ or 

1 



Then 



Ii—. 



^2 {2.557[log. 4.075B,.-1- (^+l) log. £| )' 



de 



a 



dx 



{2.667[log. 4.075B<,<,-1- (^+1 
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The variables are separable. Let 

E=2M7 log« 4.076 B,, 

J?'=2.557(g+l) 

and, when 4=1, -^=0 

The final equation then becomes 

(7|=(^^-l) (E'-2EF+2F^ 

where 

(7=1+2.036 (^^) (H-1.286) 
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F:auBE 25.— Bxponent CVC detennining rate ot piessme lecovsry for pisssme distiibntiaii 
designed to bold a constant £^ 

By the use of equations (3) and (4), curves of qjgo and 
against xJOq, which indicate the manner in which pressure 
can be recovered without a change in boundary-layer shape, 
may be plotted. Figure 26 is a plot of OijG against H", which 
shows the amoimt of pressure that can be recovered for a 
given change in as a fimction of the shape profile to be 
maintained. The plot indicates that pressure can be recov- 
ered at the most rapid rate for a value of H fairly dose to the 
value for separation, that is, about 2.3. Becaiise the flow 
with such a high value of H is apt to be unsteady, a good 
compromise between steady flow and TninimnTn increase in 
B with decrease in g would seem to be a value of Bf of 1.7 
or 1.8. 

CONCLUSIONS 

An analysis of a considerable amount of data for turbxflent 
boundary layers collected from the available literature and 
from tests in the NACA two-dimensional low-turbulence 
tunnel indicated the following conclusions: 

1. The shape of all txH-bulent-boundary-layer profiles can 
be expressed as a function of a single parameter. 

2. The variables that control the development of the 
tiurbulent boimdaiy layer apparently are (1) the ratio of the 



nondimensional pressure gradient, expressed in terms of the 
local dynamic pressme outside the boundary layer and 
boimdary-layer thickness, to the local skin-friction coefficient 
and (2) the shape of the boundary layer. 

3. An empirical equation has been developed in terms of 
these variables that, when used with the momentum equation 
and the skin-friction relation, makes it possible to trace the 
development of the turbulent boimdary layer to the separa- 
tion point. 

4. Ho systematic variation of the skin-friction coefficient 
with the diape parameter was indicated by the data. 

6. Separation occurs for values of the shape parameter 
greater than 1.8 and less than 2.6. 
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TABLE I.— CALCULATIONS OP 6 AND H FOR NACA 66(216)-222 (APPROX.) AIRFOIL SECTION 

[o-IOJ*; Ji-2.MX10«1 



lie 








c dg^ 
{• ix 


Rt 


pUi 


fi-)-2 til 
3 s<b 


i$ 
ix 




9iq79 
nixn 


24)36 
(^-1.286) 


* dr 


* ix 




0.076 


0, 674X10^ 


1.6640 


X660 


— Z490 


2.383 










0.3080 


0.666 








.086 


.6022 


1.66794 


2.537 


—2.386 


2,488 


1.800 


— LOIQ 


2.S10 


11.24 


.31B2 


.653 


—.0903115 


^517 


—.02068 


.126 


.7148 


L63726 


2.448 


—2.090 


2,903 


L74fi 


—1.082 


2,827 


2.827 


— !3606 


*5115 


—! 0001931 


— !2702 


— ! 002702 


.136 


.7429 


1.6346 


i420 


—2.066 


3,005 


1.730 


— L.I20 


2.860 


4.276 


—.8864 


!6069 


—! 0001660 


— !2234 


—! 003361 


.160 


.7867 


1.62837 


X390 


—2.040 


3,168 


L708 


— 1. 183 


2.891 


4.338 


—.3926 


.4933 


—.0001168 


— . 1474 


—.002211 


.166 


.8291 


1.62616 


2.380 


— ZOlO 


3,312 


L695 


—1.246 


ZOlO 


ia29 


—.4166 


.4888 


—.0000824 


—.09938 


— . 003478 


.200 


.9320 


1.6227 


2.290 


—1.990 


3,667 


L658 


—1.427 


8.086 


6l17 


—.4885 


.4817 


.0000077 


.00828 


.000165 


.220 


.9937 


1.62288 


2.260 


—1.940 


3,876 


L639 


— 1.C10 


3.149 


0.46 


—.5229 


.4821 


.0000465 


.0468 


.001401 


,260 


1.0882 


1.62426 


2.190 


—1.800 


4,187 


1.612 


—1.675 


3.187 


15.935 


—.665 


.4849 


.0000700 


.0734 


.003870 


.300 


i.aig 


1.62703 


2.110 


—1.600 


4,713 


L675 


— 1I68O 


3! 244 


10l220 


—1 601 


.492 


! 0001254 


!l005 


.006025 


.360 


1.410 


1.63296 


2.022 


— 1.8S0 


5,213 


LM2 


— Z316 


8.868 


3.858 


—i860 


!603 


*0QO4O0 


.2901 


.002001 


,38 


1.449 


1.63(S6 


2.005 


—2.000 


6,335 


1.636 


—2.55 


4.085 


8.170 


—.941 


.609 


.000518 


.3575 


.007160 


.38 
,40 


1.6307' 


1.64301 


1.960 


—2.240 


6, 672 


1.622 


—3.008 


4.620 


9.240 


—1.149 


.623 


.000776 


.6070 


, 01014 


1.6231 


L 66316 


1.910 


—2.530 


5,832 


l!610 


—3.820 


5.330 


10! 660 


L4a4 


I544 


.001148 


.707 


.01414 


,42 


1.7297 


1.66729 


1.88 


—2.800 


6,133 


1.495 


—4.646 


6.140 


6.140 


—1-742 


.672 


1 001615 


1.030 


.01030 


,43 


1,7911 


1.67760 


1.831 


—2.960 


6,301 


1.488 


—6.170 


0.067 


6.007 


—1.946 


l693 


.001975 


1.103 


.01103 


,44 


1,8678 


1.68862 


1.800 


—3.105 


6,481 


1.477 


—6.76 


7,227 


7,327 


—2.170 


.616 


,002393 


1.2881 


,012881 


,46 


1.9369 


LCOI60 


1.770 


—3.290 


6,700 


L47fl 


— o!488 


7! 953 


7I958 


— X449 


i642 


[002945 


1.620 


.01520 


,40 


2.0161 


1.61670 


1.731 


—3. 465 


6,898 


1.4M 


—7.279 


S.734 


8.734 


— Z766 


.673 


.003634 


1.827 


.01827 


.47 


2.104 


1.63497 


1.700 


—3.620 


7 132 


L448 


— & 142 


0.500 


4.795 


—3.091 


.710 


.004577 


Z175 


• 01088 


,476 


2.1620 


1.64686 


1.682 


—3.700 


7I267 


1.442 


—8.630 


10.075 


5.038 


3,276 


.733 


.OOSU 


2.375 


.01188 


,48 


2.2024 


1.66773 


1.663 


—3.800 


7)384 


1.441 


—0.100 


10.640 


5.320 


—3.193 


1 756 


'005S8 


2.661 


! 01330 


,486 


2.2666 


1.67103 


1.646 


—3.880 


7 622 


1.435 


—0.766 


11.200 


5.600 


—3,707 




.00661 


2.030 


.01465 


,49 


2.3116 


1.68668 


1 626 


— siOTO 


7! 661 


1.430 


—10.41 


1L840 


5,920 


—3.040 


'813 


! 007589 


3.283 


! 01642 


.406 


2.3703 


1,70210 


1.608 


-1.080 


7,816 


1.423 


-11.14 


12.56 


6.280 


—4.233 


!847 


'.00837 


3.741 


.01870 


.60 


3 4338 


1,72080 








1 AOfi 


—ll.VL 




JQ.OOO 




QQK 
. ooO 


. 01036 


^ Mil 




,606 


2.6002 


1.74208 


1.666 


-4.222 


8,136 


1.413 


-12.61 


14.03 


7.010 


-4.770 


.928 


.01210 


4.840 


.02420 


,61 


Z670 


1.76628 


1.646 


-4.262 


8,306 


1.406 


-13.36 


14.76 


7.380 


-5.043 


.977 


.01431 


5.688 


.02781 


.615 


2.6438 


1.79412 


1.626 


-4.282 


8,489 


1.401 


-14.08 


16.48 


7.740 


-6.296 


1.034 


.01709 


6.461 


.03232 


.62 


2.7212 


1.82044 


L500 


-4.202 


8,666 


L3fl6 


-14.90 


16.30 


8.160 


-6.584 


LIDO 


.02085 


7.682 


.03831 


.625 


2.8027 


1.86496 


1.476 


—4.290 


8,860 


1.391 


-16.76 


17.14. 


8.670 


-6.869 


L178 


.02621 


9.352 


.01676 


.630 


2.88S4 


1.91171 


1.450 


-4.266 


9,043 


1.384 ' 


-16,63 


18.00 


9.002 


-6.142 


1.273 


.03408 


11.7989 


.05899 
.07386 


.636 


2.9784 


1.97070 


L431 


-4.080 


9,284 


1.876 


-16.88 


18.74 


9.370 


-6.176 


1.893 


.04399 


14.770 


.64 


3,0721 


2.04466 


1.407 


-3,380 


9,474 


1.373 


-14.92 


16.29 


8.146 


-6.376 


L644 


.05019 


16.337 


.08168 


.646 


3.J638 


Z 12823 


1.390 


-3.000 


9,667 


L368 


-14.05 


16.42 


7.710 


-4.978 


1.710 


.06242 


m793 


.09896 


.66 


3.2307 


2.22610 


1.374 


-2.660 


9,846 


L362 


-13.16 


14.62 


7.280 


-4. S76 


1.911 


.07902 


24.738 


.12389 


.666 


3.3033 


2.34888 


1.362 


-2.340 


10,023 


1.367 


-12.84 


13.70 


6.850 


—4.182 


2.163 


.10681 


32.334 


.16167 


.66 
.666 


3.3718 
3.4370 


2.61066 
2.73493 


1.361 


-2.0S0 


10,201 


1.3S2 


-11.68 


13.03 


6. SIS 


-3.831 


3.192 


.16131 


14.876 


.22433 






1 



















